Fluctuation and relaxation properties of pulled fronts: 
a possible scenario for non-Kardar-Parisi-Zhang behavior 
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We argue that while fluctuating fronts propagating into an unstable state should be in the stan- 
dard KPZ universality class when they are pushed, they should not when they are pulled: The 
universal 1/t velocity relaxation of deterministic pulled fronts makes it unlikely that the KPZ equa- 
tion is the appropriate effective long-wavelength low-frequency theory in this regime. Simulations 
in 2D confirm the proposed scenario, and yield exponents (3 ~ 0.29 ± 0.01, £ ~ 0.40 ± 0.02 for 
fluctuating pulled fronts, instead of the KPZ values j3 = 1/3, £ = 1/2. Our value of (3 is consistent 
with an earlier result of Riordan et al. 



PACS numbers: 5.40+j, 5.70.Ln, 61.50.Cj 

Over a decade ago, Kardar, Parisi and Zhang (KPZ) 
I introduced their celebrated stochastic equation 
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to describe the fluctuation properties of growing inter- 
faces with height h under the influence of the noise term 
i]. A clear "derivation" of the KPZ equation is difficult 
to give, just as much as the Landau-Ginzburg- Wilson 
Hamiltonian can not straightforwardly be "derived" from 
the Ising model. However, one expects the KPZ equation 
to be the proper effective long wavelength low frequency 
theory for interfacial growth phenomena whose determin- 
istic macroscopic evolution equation is of the form 



dh 
Iff 



v(\7h) + curvature corrections. 



(3) 



Here v(Wh) is the deterministic growth velocity of a pla- 
nar interface as a function of the orientation V/i. For, 
as long as the curvature corrections of the form V 2 /i are 
nonzero, the long wavelength expansion of (^|) immedi- 
ately yields the gradient term in (Q). The philosophy is 
then that in the presence of noise, all the relevant terms 
in the KPZ equation ([[]) are generated, and that this is 
sufficient to yield the asymptotic KPZ scaling. In agree- 
ment with this picture, many interface growth models 
have been found to show the universal asymptotic 
scaling properties predicted by (|l|). 

A dynamical interface equation of the form (j|) is ap- 
propriate for interfaces whose long wavelength and slow 
time dynamics is essentially local in space and time, i.e., 
dependent on the local and instantaneous values of the 
slope and curvature. The applicability of the KPZ equa- 
tion is therefore not limited to situations with a micro- 
scopically sharp interface: Many pattern forming systems 
of the reaction- diffusion type exhibit fronts whose intrin- 
sic width / is finite. For curvatures k small compared to 



I ~ 1 , kI <C 1 , an effective interface approximation or mov- 
ing boundary approximation of the form ^) can then be 
derived using standard techniques || . These approxima- 
tions apply whenever the internal stability modes of the 
fronts relax exponentially on a short time scale, so that an 
adiabatic decoupling becomes exact in the limit nl — > 0. 
The best known example of such a type of analysis is for 
the curvature driven growth in the Cahn-Hilliard equa- 
tion, but moving boundary techniques have recently been 
applied successfully to many other such problems |^|. In 
all these cases, the internal relaxation modes within the 
fronts or transition zones are indeed exponentially decay- 
ing on a short time scale. 

From the above perspective, recent results for the re- 
laxation properties of planar fronts propagating into an 
unstable state, suggest an interesting new scenario for 
non-KPZ behavior. Fronts propagating into unstable 
states generally come in two classes, so-called pushed 
fronts and pulled fronts 0. Pushed fronts propagat- 
ing into an unstable state are the immediate analogue 
of fronts between two linearly stable states. In the thin 
interface limit, kI <C 1, the dynamics of such fronts 
becomes essentially local and instantaneous, and given 
by an equation of the form (|^); according to the argu- 
ments given above, fluctuating pushed fronts should thus 
obey KPZ scaling: following standard practice by say- 
ing that the d+lD KPZ equation (where the +1 refers 
to the time dimension) describes the fluctuations of a d- 
dimensional interface, the conclusion is that fluctuations 
of d-dimensional pushed fronts in (d+1) bulk dimensions 
are described by the d+lD KPZ equation. 

Pulled fronts, however, behave very differently from 
pushed ones. A pulled front propagating into a linearly 
unstable state is basically "pulled along" by the linear 
growth dynamics of small perturbations spreading into 
the linearly unstable state. The crucial new insight for 
our discussion is the recent finding that pulled fronts 
can not be described by an effective interface equation like 
(0) that is local and instantaneous in space and time, even 
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if they are weakly curved on the spatial scale. This just 
reflects the fact that the dynamically important region 
of pulled fronts is the semi-infinite leading edge region 
ahead of the front, not the nonlinear front region itself. 
Technically, the breakdown of an intcrfacial description is 
seen from the divergence of the solvability type integrals 
that arise in the derivation of a moving boundary approx- 
imation in dimensions d>2 More intuitively, the re- 
sult can be understood as follows: a deterministic pulled 
front in d=l relaxes to its asymptotic speed v* with a uni- 
versal power law as v (t) = v*+ Cl /t+c 3/2 /t 3 / 2 + - ■ where 
ci(< 0) and c 3 / 2 are known coefficients ||. Clearly, this 
very slow power law relaxation implies that an adiabatic 
decoupling of the internal front dynamics and the large 
scale pattern dynamics can not be made and hence that 
there is no long- wavelength effective interface equation of 
the form (||) for pulled fronts. There is then a priori no 
reason to expect that fluctuating pulled fronts are in the 
KPZ universality class! 

It is our aim to test this scenario by introducing a sim- 
ple stochastic lattice model whose front dynamics can be 
changed from pushed to pulled by tuning a single param- 
eter. Our results are consistent with our conjecture that 
pulled fluctuating fronts are not in the standard KPZ 
universality class, while pushed fronts are. In fact, our 
results put an earlier empirical finding of Riordan et al. 
H into a new perspective: These authors obtained es- 
sentially the same growth exponent as we do for the non- 
KPZ case, but the connection with the transition from 
pushed to pulled front dynamics was not made. 



for all values k% < yk\k$/2 and pushed for larger k 2 . 
In the pulled regime, the asymptotic front velocity, is 
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(a) Pushed 



(b) Pulled 



FIG. 1. The two versions of our stochastic model for noisy 
pushed and pulled fronts. The model is isotropic, i.e., all 
neighbors are probed with equal probability, (a) The pushed 
case: Stochastic moves consist of diffusive jumps of a particle 
to a neighboring empty site and birth and death processes at 
sites whose two randomly chosen nearest neighbors are oc- 
cupied, (b) The pulled case: The only difference with the 
pushed case is in the birth and death processes. 

Our stochastic model is motivated by || and the re- 
sults for deterministic planar fronts in the nonlinear dif- 
fusion equation 



dp/dt = DW 2 p + hp + k 2 p 2 - k 3 p 3 



(4) 



As discussed in 10 7|], the planar fronts with p > 
propagating into the unstable state p = are pulled 



v* = 2y/Dki, while in the pushed regime the asymptotic 
front velocity equals 
where K 



2VDh[(-K + VA' 2 + 4)/32] 
We confine ourselves here to 



studying two limits where the stochastic front dynam- 
ics can easily be understood intuitively. 

We study the dynamics of particles on a square lattice, 
subject to the constraints that no more than one parti- 
cle can occupy each lattice site. The stochastic moves are 
illustrated in Fig. 1. They consist of diffusive hops of par- 
ticles to neighboring empty sites and of birth and death 
processes on sites neighboring an occupied site. In a 
mean field approximation, this stochastic model is equiv- 
alent to a discrete version of (^) . We will study here the 
two cases indicated in Fig. 1. For k\ — (Fig. la), planar 
fronts are definitely pushed: Since the linear spreading 
speed v* = for k\ = 0, the front must then be pushed, 
even if corrections to the mean field behavior are impor- 
tant in the front region or behind the front. Likewise, 
when k 3 = and k 2 < (Fig. lb), the nonlinearities 
behind the front only limit the birth (growth) rate, so in 
this limit the stochastic planar front is definitely pulled. 




i=L 



FIG. 2. A snapshot of the coarse-grained density field 
(m = 1). The interface position (continuous line) is obtained 
as the point where the density crosses half its bulk value. 

Our simulations are done on 2D strips which are long 
in the y direction and of width L in the x direction. In 
the x direction, periodic boundary conditions are used. 
The Monte Carlo simulations are started with a configu- 
ration in which the first few rows (~ 100) of the lattice 
are occupied with a probability equal to the equilibrium 
density. All other lattice sites are empty. After an initial 
transient, the scaling properties of the interface width are 
studied in the standard way using the following defini- 
tion of the interface height h. We define a coarse-grained 
density variable at each lattice site as the average occu- 
pation of sites on a (2m + 1) x (2m + 1) grid centered at 
that site. We then define the position h(xi) of the inter- 
face as the first point where this coarse-grained density 
reaches half the equilibrium density value. Our results 
for ensemble averaged width of the interface (see below) 
are obtained by averaging over 100 runs for the largest 
system L = 2048 to about 3200 runs for the smallest 
L = 64. Although we have performed simulations with 
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m = 1,2 and m = 3 almost all the data presented sub- 
sequently are those for a representative value of m = 2. 
The coarse-grained density field and the corresponding 
interface position h for a typical configuration is shown 
in Fig. 2. 
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FIG. 3. Scaling collapse of the width in the pushed (empty 
symbols) and the pulled (solid symbols) cases using (a) KPZ 
exponents ( = 0.5,2 = 1.5 and (b) £ = 0.4, z = 1.38. The 
pushed data values are multiplied by a factor 3 for clarity. 

The interface width it; of a given realization is defined 

in the usual way, w 2 (t) = (h(xi,t) — h(xi,t)) 2 , where the 
overbar denotes a spatial average, h = L^ 1 ^2 X . h(xi,t). 
The proper scaling to study is the ensemble averaged 
mean square interface width W 2 = (w 2 ). As is well 
known, in the KPZ equation W obeys a scaling form 
W(t) = tPy(j?)- Here the scaling function y(u) is 
about constant for u <C 1 and y ~ u^ 13 for u 3> 1, 
with the KPZ exponents z = 3/2 and j3 = 1/3 in 1 + ID. 
For t ^S> L z , the width saturates at W sa t ~ where 
£ = (3z is the roughness exponent. In Fig. 3, we show our 
data for stochastic pushed and pulled fronts by plotting 
W/L^ versus t/L z for a range of system sizes (L — 128 to 
2048). Following standard practice, we always plot the 
subtracted width W 2 (t) — !F 2 (0) to minimize the effect 
of the initial front width. The kinetic parameters are 



chosen to be k% = 0.5, k% = 1.0 for the pushed model and 
fc 2 = 0, ki =0.1 for the pulled model. The diffusion rate 
D = 0.25 is the same in both cases. In Fig. 3a we use the 
1+1D KPZ exponents to obtain a data collapse. Clearly, 
good scaling collapse of the pushed data confirms that the 
pushed fronts are in the universality class of the 1+1D 
KPZ equation. By contrast, use of 1+1D KPZ values 
does not lead to good scaling collapse of the pulled data. 
In Fig. 3b we show the same sets of data but now with 
exponents £ = 0.4 and z = 1.38 to obtain the best possi- 
ble scaling collapse of the pulled data. It is clear that the 
two sets of exponents, though only moderately different 
from each other, are well beyond error-bars. More accu- 
rate estimates of the exponents for the pulled case were 
obtained as follows. In Fig. 4 we fit a power law to the 
non-saturated part of the width for the largest system 
L = 2048 and obtain j3 ~ 0.29 ± 0.01 for the growth ex- 
ponent. Plotting the saturated width W sa t as a function 
of system size L (Fig. 4, inset) yields ( ~ 0.4 ± 0.02 for 
the roughness exponent. Once £ is known the dynamic 
exponent is obtained by requiring good scaling collapse 
of Fig. 3b., z — 1.38 ± 0.06. The value of /3 is consistent 
with that reported by Riordan et al. [B for this model, 
(3 = 0.272 ± 0.007, but their apparent value of z ~ 1 
is not the true dynamic exponent related to the inter- 
face roughness through £ = j3z, since they studied the 
ensemble averaged width of the front [|ll] . 
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FIG. 4. Scaling of the width in the pulled case for a system 
of size L — 2048. We have taken ki = 0.1, /c2 = —1.0 and 
D — 0.25. The inset shows the saturated width W sa t vs L 
plot on a log-log scale for L — 64 x 2 n (n = 0, ■ ■ ■ , 5) and is 
consistent with £ = 0.4 (solid line). 

Another way to investigate the possible difference with 
the 1+1D KPZ behavior is to study the distribution 
P(w 2 /W 2 ). For ID interface models whose long time 
interface configurations are given by a Gaussian distri- 
bution, like the KPZ model, the distribution function 
P(w 2 /W 2 ) is uniquely determined, without adjustable 
parameters Jl3[ . As Fig. 5 shows, in the pushed regime 
our data are completely consistent with this distribution 
function, but in the pulled regime the measured distri- 
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bution function deviates significantly from the universal 
prediction for Gaussian interface fluctuations. 
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FIG. 5. The probability distribution of the width of the 
interface (a; = w 2 /W 2 ). The triangles and the circles are for 
the pushed and the pulled data of Fig. 3, respectively, and 
the solid and the dashed lines are the universal distribution 
functions for the 1 + 1D KPZ and 2 + ID KPZ equation, as ob- 
tained by Racz et al. |l3"| . Inset: the same plot on a semi- log 
scale showing the agreement for large arguments. 

The essential difference between pushed and pulled 
fronts is that for pushed fronts the dynamically impor- 
tant region is the finite transition zone between the two 
phases it separates, whereas for pulled fronts it is the 
semi-infinite leading edge ahead of the front itself |7],|8| . It 
is precisely for this reason that the wandering of stochas- 
tic pulled fronts in one bulk dimension with multiplica- 
tive noise was recently found to be subdiffusive and de- 
termined by the 1+1D KPZ equation, not by a "0+1D" 
stochastic Langevin equation [|12| . By extending this idea 
it has been recently conjectured [B that the scaling ex- 
ponents of stochastic pulled fronts in d+1 bulk dimen- 
sions are generally given by the (d+l)+lD KPZ equation 
instead of the d+lD KPZ equation, essentially because 
the dimension perpendicular to the front can not be inte- 
grated out [Q. The scaling exponents we find here in 2 
bulk dimensions are indeed close to those reported for the 
2+1D KPZ equation M, the supposedly most accurate 
values being C = 0.393(3), f3 = 0.245(3) ||. Moreover, 
the probability distribution P{w 2 /W 2 ) of pulled fronts 
fits the P{w 2 /W 2 ) of the 2+1D KPZ equation quite well 
without adjustable parameters, see Fig. [|. For justifica- 
tion and further exploration of this conjecture, we refer 
to @. 

An interesting limit of our model is obtained when we 
further take k% — in Fig. lb. In this case only birth and 
diffusion occurs, leading to an equilibrium density p = 1 
behind the front. If we put D = as well, the result 
is an Eden-like model j4|^] with the modification that 
the probability of adding a particle is proportional to the 



number of neighbors, not independent of it. Numerical 
simulations in this Eden-like limit indicate that the KPZ 
exponents are recovered, as it should, and hence that the 
model has a KPZ to non-KPZ transition at intermediate 
values of p eq and D. 

In conclusion, even though one should always be aware 
of the possibility of a very slow crossover to asymp- 
totic behavior in such studies |l6| ] - a problem that 
has plagued some earlier tests of KPZ scaling in, e.g., 
the Eden model — taken together our data as well as 
those of H give, in our opinion, reasonably convincing 
evidence for our scenario that the absence of an effective 
interface description for deterministic pulled fronts also 
entails non-KPZ scaling of stochastic pulled fronts. 
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